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Unidirectional and robust transport is generally observed at the edge of two- or three-dimensional 
quantum Hall and topological insulator systems. A hallmark of these systems is topological pro¬ 
tection, i.e. the existence of propagative edge states that cannot be scattered by imperfections or 
disorder in the system. A different and less explored form of robust transport arises in non-Hermitian 
systems in the presence of an imaginary gauge field. As compared to topologically-protected trans¬ 
port in quantum Hall and topological insulator systems, robust non-Hermitian transport can be 
observed in lower dimensional (i.e. one dimensional) systems. In this work the transport properties 
of one-dimensional tight-binding lattices with an imaginary gauge held are theoretically investigated, 
and the physical mechanism underlying robust one-way transport is highlighted. Back scattering 
is here forbidden because reflected waves are evanescent rather than propagative. Remarkably, the 
spectral transmission of the non-Hermitian lattice is shown to be mapped into the one of the cor¬ 
responding Hermitian lattice, i.e. without the gauge held, but computed in the complex plane. In 
particular, at large values of the gauge held the spectral transmittance becomes equal to one, even 
in the presence of disorder or lattice imperfections. This phenomenon can be referred to as one-way 
non-Hermitian transparency. Robust one-way transport can be also realized in a more realistic set¬ 
ting, namely in heterostructure systems, in which a non-Hermitian disordered lattice is embedded 
between two homogeneous Hermitian lattices. Such a double heterostructure realizes asymmet¬ 
ric (non-reciprocal) wave transmission. A physical implementation of non-Hermtian transparency, 
based on light transport in a chain of optical microring resonators, is suggested. 

PACS numbers: 03.65.-w, 72.20.Ee, 72.15.Rn, 73.43.-f, 71.30.+h 


I. INTRODUCTION 

Electrons in two-dimensional structures with a mag¬ 
netic field exhibit a wide variety of collective quantum 
phenomena such as integer, fractional and quantum spin 
Hall effects CHS- A hallmark of these systems is the 
existence of edge states with topological protection that 
can propagate at the boundaries immune of disorder and 
scattering. The robustness of topological states has been 
actively explored in a wide variety of quantum Hall sys¬ 
tems, topological insulators, and topological supercon¬ 
ductors (see e.g. [5] [6] and references therein), with 
potential applications in metrology, quantum computing 
and spintronics. Recently, a large amount of theoreti¬ 
cal and experimental effort has been devoted toward the 
emulation of quantum Hall behavior in bosonic systems, 
i.e. in ultra-cold gases mM an( l photons [SHU]. Be¬ 
sides of the possibility to explore new quantum phases 
of matter, in optics topological insulation and protec¬ 
tion hold the promise for applications in optical isola¬ 
tion and robust photon transport. However, propagative 
states with topological protection do not arise in one¬ 
dimensional (ID) systems. A natural question thus arises 
whether quantum-Hall-like robust transport, which is in¬ 
sensitive to disorder and imperfections, can be realized in 
ID systems. In 1996, Hatano and Nelson [2T] investigated 
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the problem of Anderson localization in a ID disordered 
non-Hermitian lattice. They showed that an ’imaginary’ 
(rather than real) magnetic field can prevent Anderson 
localization, with the appearance of a mobility interval at 
the center of the band. Such a result, referred to as non- 
Hermitian delocalization, has been subsequently revisited 
by several authors [22143T] and raised some debate about 
the nature of eigenstates. Inspired by the Hatano-Nelson 
delocalization transition, in a recent work [521 it has been 
conjectured that robust transport is expected to arise in 
a rather general class of non-Hermitian ID lattices, pro¬ 
vided that a simple condition for the lattice energy band 
is satisfied. The proposal of an optical implementation 
of an ’imaginary’ gauge field, based on a chain of coupled 
optical microrings with tailored gain and loss regions I32j , 
renewed the interest in the Hatano-Nelson model. 

In this paper we present a theoretical study of robust 
one-way transport in ID tight-binding lattices with an 
imaginary gauge field, highlighting the phenomenon of 
non-Hermitian transparency: besides to prevent Ander¬ 
son localization and restoring mobility, it is shown that 
in the presence of a sufficiently high imaginary magnetic 
field structural imperfections and disorder in the lattice 
become almost one-way transparent , i.e. they do not scat¬ 
ter back waves and do not substantially alter the spectral 
transmission of the lattice. Such a result follows from the 
fact that (i) back-scattered waves are evanescent waves, 
and (ii) the spectral transmission of the non-Hermitian 
lattice can be mapped into the one of the correspond- 
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ing Hermitian lattice, i.e. without the imaginary gauge 
field, but computed in the complex energy plane. In par¬ 
ticular, at large values of the gauge field the spectral 
transmittance becomes equal to one, even in the presence 
of disorder or lattice imperfections, thus realizing non- 
Hermitian transparency. Robust one-way transport and 
transparency effects can be also realized when the imag¬ 
inary magnetic field is applied to a limited spatial region 
of the array. In this case, a suitable double heterostruc¬ 
ture system, in which a non-Hermitian disordered lattice 
is embedded between two homogeneous Hermitian lat¬ 
tices, can realize one-way transparency. 


II. SCATTERING IN A TIGHT-BINDING 
LATTICE WITH AN IMAGINARY GAUGE 
FIELD: ROBUST TRANSPORT AND ONE-WAY 
NON-HERMITIAN TRANSPARENCY 


A. Lattice model 
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We consider the hopping motion of a quantum particle 
in a disordered ID tight-binding lattice in the presence 
of an imaginary gauge field [2 Tj, see Fig. 1(a). In the 
nearest-neighbor approximation, the system is described 
rather generally by the non-Hermitian Hamiltonian 

H = ^2 ( K ™-i,n\ n - l)(n| + K n ,n-i \n)(n - 1|) 

n 

+ ^2 V n\n)(n\ ( 1 ) 


where |n) is a Wannier state localized at site n of the 
lattice, K n , n - t-i and K n +i ;n are the generally asymmetric 
(i.e Kn t n +1 7 ^ k* + 1 „) hopping rates between site | n) and 
|n + 1), and V n is the energy of Wannier state |n). After 
setting = J2n c n(t)\ n ) an d assuming h = 1, the 

amplitude probabilities c n {t) satisfy the set of coupled 
equations 


,dc r 

l-r- 
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( 2 ) 


where we have set 


Gn,m "b ^n,n— l^n— l,m T V n 5 n , r 


(3) 


The energy spectrum E of TL is obtained from the eigen¬ 
values of the tridiagonal matrix Q, i.e. as eigenvalues of 
the linear equations 


E'lpn — ^n.n+lt^n+l d" ^n 7 n—l' i Pn — l T V n 1p n (4) 

with appropriate boundary conditions. In the Hatano- 
Nelson model [2T], disorder is allowed for the site energies 
V n , whereas the hopping rates are given by 

K n ,n+ 1 = ft exp (~h) , = «exp(/i) (5) 


FIG. 1. (color online) (a) Schematic of a one-dimensional 
non-Hermitian tight-binding lattice with asymmetric hop¬ 
ping rates induced by an imaginary gauge field h [«„,„+1 = 
/texp(— h), Kn,n-i = Kexp(h)]. (b) Optical realization of the 
non-Hermitian lattice model shown in (a) based on photon 
hopping in coupled microring resonators, (c) Schematic of 
a double heterojunction lattice: a disordered lattice region, 
with site energy disorder W„ and applied imaginary gauge 
field h, is connected at sites n = —N and n = N to two 
Hermitian homogeneous lattices. 


where h > 0 describes the effect of the imaginary vec¬ 
tor potential and k is the hopping rate in the absence 
of the imaginary gauge field. The Hermitian lattice case 
is obtained in the limit h 0. The possibility to im¬ 
plement asymmetric hopping rates as given by Eq.(5) 
was originally discussed by Hatano and Nelson for mag¬ 
netic flux lines in type-II superconductors [25]. In op¬ 
tics, realization of asymmetric hopping rates was sug¬ 
gested in a few different settings, including engineered 
coupled waveguide lattices [55!, active mode-locked lasers 
with combined amplitude and phase modulators [33| , and 
chains of coupled optical microring resonators [32]. Fig¬ 
ure 1(b) shows a possible experimental implementation 
of an imaginary gauge field ft in a chain of coupled mi¬ 
croring optical cavities, which has been recently proposed 
in Ref. [32j. The resonator chain consists of a sequence of 
main ring resonators which are indirectly coupled using 
an interleaved set of auxiliary rings. The auxiliary rings 
are designed to be antiresonant to the main ring res¬ 
onators, i.e., the length of the connecting rings is slightly 
larger (or smaller) than the main rings so as to acquire 
an extra 7r phase shift. To realize a synthetic imaginary 
gauge field h, the auxiliary ring provides amplification 
in the upper half perimeter, with single-pass amplifica¬ 
tion h , and balanced loss in the lower half perimeter, with 
single-pass attenuation —h. Light circulation in the main 
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(auxiliary) rings is forced to be counterclockwise (clock¬ 
wise) by some non-reciprocal element or just because of 
the excitation conditions. Indicating by c n (t) the am¬ 
plitude of the counterclockwise propagating field in the 
n-th ring in the main resonators, with a carrier frequency 
coincident with one longitudinal ring resonance, coupled¬ 
mode equations for the slowly-varying amplitudes c n (t) 
can be derived in the mean-field limit and are precisely 
of the form (2), with hopping rates given by Eq.(5) 152] , 
The explicit expression of the hopping rate k in terms 
of coupling constants of the microrings can be found in 
Ref. [T?| . In the optical microring structure of Fig. 1(b) 
the disorder V n arises because of deviations of microring 
resonance frequencies from the reference value due to fab¬ 
rication imperfections, as well as from slight deviations 
of the antiresonance condition in the auxiliary rings ?>2| . 

To properly define the scattering process induced by 
disorder, we consider an infinitely extended lattice and 
assume that disorder or lattice imperfections are confined 
in the region —N <n<Noi lattice sites, i.e. we assume 

V n = 0 for \n\ > N. (6) 

Let us first consider the case of a homogeneous lattice, 
i.e. V n = 0. In this case the (improper) Bloch-Floquet 
eigenfunctions of H are given by il> n = exp (—iqn) where q 
is the Bloch quasi momentum that varies in the first Bril- 
louin zone (—n < q < n). The corresponding eigenener- 
gies E form a tight-binding band with complex energy 
spectrum described by the dispersion curve 

E(q) = 2 k cos(q — ih) (7) 

= 2k cosh (h) cos (q) + 2*Ksinh(/i) sin(g). 

The group velocity of a wave packet (also related to the 
current EH 1231) , obtained as a superposition of Bloch- 
Floquet states with carrier wave number q , is given by 
v g = —R e(dE/dq) = — 2«cosh(/i) sin(g). Note that a 
forward-propagating wave (0 < q < 7r, v g > 0) is ampli¬ 
fied because Im(i?(g)) > 0, whereas a backward propa¬ 
gating wave (—7r < q < 0, v g < 0) is attenuated because 
lm(E(q)) < 0. As discussed in Ref.E2, this circumstance 
makes wave transport in the lattice highly asymmetric 
because backward-propagating waves vanish after some 
propagation distance due to damping whereas forward- 
propagating waves are amplified. A different and may 
be simpler point of view to understand asymmetric (one¬ 
way) transport in the lattice when ft / 0 is the following 
one. The imaginary gauge field h increases the effective 
hopping rate in one direction by a factor exp(/i), while 
in the opposite direction it is diminished by the factor 
exp(—fti), i.e. 2 h represents the exponential suppression 
of the rates in the two directions. Such an asymmetry 
in left/right hopping rates leads to a preferred trans¬ 
port in the direction of increased hopping rates, and it 
is thus similar to asymmetric transport investigated in 
other physical contexts (see, for example, |34)[35j). Note 
that, by reversing the sign of h 7 the direction of trans¬ 
parent propagation is reversed. 


(a) (b) 



lattice site n lattice site n 


FIG. 2. (color online) Propagation of a Gaussian wave packet 
(snapshots of |cn(t)| 2 /l c ™W| 2 ) h 1 a uniform lattice with 
site-energy disorder for (a) h = 0 (Hermitian lattice), and 
(b) h = 0.2. The energies V n at sites —N < n < N (with 
N = 10) are taken from a uniform distribution in the range 
(—2k, 2k) (see the inset). Upper panels refer to left-side 
incidence (carrier Bloch wave number of the incident wave 
packet q = 7r/2), whereas lower panels refer to right-side inci¬ 
dence (carrier Bloch wave number of the incident wave packet 
q = -7t/2). 


B. Scattering states and Non-Hermitian 
transparency 

Extended numerical simulations of wave packet prop¬ 
agation in the Hatano-Nelson lattice showed that trans¬ 
port is robust against disorder and lattice imperfections 
for forward-propagating waves, but not for backward 
propagating waves l32| . An example of this behavior is 
shown in Fig.2. Such a result was attributed to the asym¬ 
metry of the complex energy dispersion curve for for¬ 
ward and backward propagating waves, and was related 
to the non-Hermitian delocalization transition predicted 
in the original work by Hatano and Nelson. However, 
the general properties of scattering states of the lattice 
and the spectral transmittance for forward waves in the 
presence of disorder were not investigated in such previ¬ 
ous works. Here we close this gap and provide analytical 
form for the spectral transmission function of forward- 
propagating waves, predicting an important phenomenon 
that was not disclosed in previous works: non-Hermitian 
transparency. This means that, for a sufficiently large 
imaginary gauge field, besides to prevent Anderson lo¬ 
calization the transport properties in the lattice are com¬ 
pletely insensitive to the disorder, i.e. the spectral trans¬ 
mission, both in phase and amplitude, is not altered by 
the presence of disorder or lattice imperfections. 

Let us consider a lattice with Hermitian disorder or im¬ 
perfections of lattice sites V n confined in the region — N < 
n < N [see Eq.(6)]. We look for a scattered state solu¬ 
tion %l) to Eq.(4) with energy E = E(q) = 2kcos(<7 — ih) 
(0 < q < 7r), corresponding to a forward propagating 
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FIG. 3. (color online) Numerically-computed spectral transmission t(q) (amplitude square and phase) for left-side incidence 
in a uniform lattice with site-energy disorder at sites — N < n < N with N = 20 for (a) h = 0 (Hermitian lattice) and (b) 
for a few increasing values of the gauge field h (curve 1: h = 3, curve 2: h = 5, curve 3: h = 7). The energies V„ are taken 
from a uniform distribution in the range (—2k, 2k) [see the inset in (a)]. Panel (c) shows the numerically-computed spectral 
transmission (amplitude square and phase) for h = 5 and for 20 different realizations of lattice disorder. 


wave ~ exp(— iqn) with Bloch wave number q coming 
from n = — oo and incident upon the disordered region 
of the lattice. To this aim, let us first notice that in the 
limiting case of an Hermitian lattice (h = 0) imperfec¬ 
tions and disorder generate a backward-propagating wave 
with opposite Bloch wave number q , and the asymptotic 
form of the scattered state is given by 

exp (—iqn) + r 0 (q) exp (iqn) (n < —N) , , 

to(q) exp(— iqn) ( n>N ) ' ' 


where rq(q) and to(q) are the spectral reflection and 
transmission coefficients, respectively. The explicit ex¬ 
pression of the spectral transmission and reflection coef¬ 
ficients can be determined by standard transfer matrix 
method [3B] and read (see Appendix A for technical de¬ 
tails) 


to(q) 

r o (?) 


2 i sin(g) 


Vn exp (iq) - V 22 exp (—iq) + V 12 - V 21 
Q 22 ~ Q 11 + Q 21 exp (-iq) - 0,12 exp (iq) 
Q 11 exp (iq) - Q 22 exp (-iq) + Q 12 - Q 2 1 


exp[iq(2N + 1)] 


exp[iq(2N + 1)] 


(9) 

( 10 ) 


where 

we have set 



Q(q) = M.n x Mn-i x •■■■ x ■M-n 

( 11 ) 


V{q) = M.~n x A4 _at+i x .... x A4 jv 

( 12 ) 

and 

-M„M = ( 2cos(,) y’ 4/, ‘ - 1 ). 

(13) 


conservation a scattering process of the incident wave 
with Bloch wave number q generates a wave with Bloch 
number q\ satisfying the condition E(qi) = E(q). Using 
Eq.(7), such a condition reads explicitly 

qi=-q + 2ih (14) 


In the presence of disorder and for sufficiently large values 
of N, one has to(q) —> 0 and r 0 (q) —► 1 owing to Anderson 
localization; see for example Fig.3(a). 

In the presence of the imaginary gauge field (/i ^ 0) a 
fully different scenario is found. In fact, owing to energy 


i.e. the scattered wave is an evanescent wave because it 
corresponds to a complex wave number qi. Such a result 
explains why the imaginary gauge field suppresses back 
reflections in the lattice. For h 7 ^ 0, the scattered state 
solution to Eq.(4), corresponding to left-side incidence, 
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can be searched in the form [compare with Eq.(8)] 

! exp(— iqn) + r(q) exp[2 h(n + AT)] exp (iqn) 

(n < -N) (15) 

t(q) exp {—iqn) (n > N) 

where t(q), r(q) are the spectral transmission and re¬ 
flection coefficients. A simple relation can be established 
between t(q), r(q) and their limits to(q), r 0 (q) in the Her- 
mitian case h = 0. In fact, after setting ip n = (f> n exp (hn) 
and p = q — ih , from Eqs.(4), (5) and (15) one obtains 

E(j) n = K(<j) n+1 + <j> n -i) + V n <j> n (16) 

with E = 2k cos(p) and with the asymptotic behavior 

{ exp(— ipn) + r(q) exp(2 hN) exp (—ipn) 

{n < -N ) (17) 

t(q) exp (—ipn) (n > N) 

Note that the scattering problem defined by Eqs.(16) and 
(17) is equivalent to the one of the Hermitian lattice, but 
with the real Bloch wave number q replaced by p = q—ih. 
Such an equivalence also follows from an inspection of 
the band dispersion relation (7), in which we realize that 
the lattice band of the non-Hermitian system (i.e. for a 
non-vanishing imaginary gauge field h 0) is obtained 
from the Hermitian one (i.e. in the limit h = 0) after 
complexification of the Bloch wave number according to 
the relation p = q — ih. A comparison of Eqs.(8) and (17) 
then yields 

t(q) — to(q — ih) (18) 

r(q) = exp(—2Nh)ro(q — ih). (19) 

It is worth considering the behavior of t(q) and r(q) in the 
limit of a large imaginary gauge field, i.e. for h —» oo. 
Using the expressions (9) and (10) for t 0 (q), ro(q) and 
from Eqs.(18) and (19) it can be proven that the following 
asymptotic relations hold (see Appendix B) 

t(q) 1 , r(q) —» 0 as h — > oo. (20) 

Equation (20) shows that, for a high enough imaginary 
gauge field h, forward wave packet propagation is not af¬ 
fected by disorder or imperfections in the lattice, i.e. uni¬ 
directional non-Hermitian transparency is realized. This 
phenomenon is illustrated in Figs.3(b) and (c). Figure 
3(b) shows the behavior of the spectral transmission t(q) 
(amplitude and phase) in the same disordered lattice 
for increasing values of h. Note that, according to the 
asymptotic analysis, as the magnetic flux h is increased 
the lattice becomes transparent and the spectral trans¬ 
mission t gets close to one. Figure 3(c) shows the behav¬ 
ior of the spectral transmission for a fixed value of the 
magnetic field (h = 5) and for 20 different realizations 
of disorder. The figure clearly indicates that the trans¬ 
mission is almost insensitive to the disorder realization 
in the lattice. It should be emphasized that the appear¬ 
ance of non-Hermitian transparency is closely related to 


the complexification of the Bloch wave number. In fact, 
Eq.(18) basically states that the transmission coefficient 
of propagative waves in the non-Hermitian case is ob¬ 
tained from the expression of the transmission coefficient 
t(jp) in the Hermitian case after the complexification of 
the Bloch wave number p = q — ih. Why is it important 
such a complexification to realize non-Hermitian trans¬ 
parency? The answer is as follows. For a real p, i.e. for 
h = 0, t(q) is the transmission coefficient of the Hermitian 
lattice, which is strongly sensitive to the site energy po¬ 
tentials V n : backscattering is fully in action, leading e.g. 
to Anderson localization in the disordered lattice. When 
we extend t(p) into the complex plane, t(p) becomes less 
sensitive to the site potentials V n , and full transparency 
t(p) ~ 1 is realized in the limit h —> oo, as briefly sketched 
in Appendix B. Physically, the lack of sensitivity to V n for 
a large imaginary gauge field h relies on the fact that the 
back-scattered waves are evanescent ones and the cou¬ 
pling to them is weaker and weaker as h. is increased: ba¬ 
sically in the forward-propagation direction a wave does 
not see V n as scattering centers. Such a result holds not 
only for a random distribution of site energies V n , but also 
for structural imperfections or defects in the lattice. For 
example, when V n is non-vanishing in two sites of the lat¬ 
tice, in the Hermitian lattice (h = 0) a propagating wave 
packet undergoes multiple reflections back and forth be¬ 
tween the two defects, like in a Fabry-Perot cavity. This 
yields multiple transmitted and reflected wave packets, 
i.e. echoes of the original wave packet are observed (see 
for example Fig. 3(a) in Ref. [32]). Application of the 
imaginary gauge field to the lattice (h ^ 0) makes the 
two defects invisible and echo effects are suppressed for 
forward-propagating waves. 

A fully opposite behavior is found when considering a 
wave that propagates in the backward direction. As 
shown in Fig.2(b), here strong back scattering induced 
by disorder arises, with the reflected wave packet being 
amplified. For the ideal case of an imaginary gauge field 
uniformly applied all along the entire lattice, the scatter¬ 
ing problem for right-side incidence is formally ill posed. 
The reason thereof is that the scattered states for right- 
side incidence contain reflected evanescent waves which 
are undamped as n +oo. However, such unbounded 
states can be superimposed to describe propagation of 
localized wave packets, which have physical meaning. In 
such an analysis, the unbounded evanescent waves lead 
to an effective amplification of the back scattered wave 
packet, as observed in Fig.2(b). A detailed analysis of 
such a case is not given here. Indeed, in the more realistic 
case considered in the next section, with the imaginary 
gauge field applied to a limited spatial region of the lat¬ 
tice, the scattering problem for plane waves is well posed 
for both incidence sides. 
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(a) 
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FIG. 4. (color online) (a) Numerically-computed spectral 
transmittance ( \t^ I,r \q)\ 2 ) and reflectance ( \r^ l,r \q)\ 2 ) of a 
heterostructure lattice (N = 10), for left and right sides of in¬ 
cidence, in the absence of disorder. The different curves refer 
to increasing values of the gauge field h. Curve 1: h — 0.2; 
curve 2: h = 1; curve 3: h = 1.5; curve 4: h = 2.5. The 
dashed curves show the behaviors of spectral transmittance 
and reflectance in the h —> oo limit [Eqs. (34-37) given in the 
text]. (b) Behavior of the phase of the spectral transmittance 
t«(g). 


III. DOUBLE HETEROSTRUCTURE LATTICE 
AND ONE-WAY ROBUST TRANSMISSION 

In the previous section we assumed an infinitely ex¬ 
tended lattice with a uniform imaginary gauge field h 
extended over the entire lattice, both in the ordered and 
disordered regions. However, in a more realistic system 
one wishes to apply the imaginary gauge field to a lim¬ 
ited spatial region of the lattice, namely to the one with 
disorder, with the aim to prevent Anderson localization 
and to suppress the impact of disorder on the spectral 
transmission. In this case one basically realizes a double 
heterostructure: a non-Hermitian lattice with disorder 
embedded between two Hermitian homogeneous lattices; 
see Fig. 1(c). The Hamiltonian of the heterostructure is 
described by Eq.(l) with 

K n ,n+1 = Hn+l,n = « U < -N - 1 , 71 > N 

K„,n +1 = Kexp(— h) — N < n < N — 1 (21) 

Kn } n -1 = Kexp (h) — N + 1 < n < N. (22) 

The site energy V n is assumed uniform and equal to zero 
in the Hermitian regions (V n = 0 for n < —N and n > 
N), whereas in the non-Hermitian region we now assume 

V n = -iS + W n (-N <n < N). (23) 

In Eq.(23) S provides an imaginary bias to the energy 
whereas the real parameters W n account for disorder and 


lattice imperfections in the site energies. The Hermitian 
and non-Hermitian sections are attached at the sites n = 
—N and n = N, at which the site energies are V_n an d 
V)v will be determined later. The spectral transmission 
and reflection of the structure, for both left (l) and right 
(r) incidence sides are defined as the scattering states 
[Eq.(4)] with the asymptotic behavior 


V’n = 


exp (—iqn) + A l \q) exp (iqn) (n < — N ) 


t^(q) exp (—iqn) 
for left-side incidence, and 


(n > N) 


4>n = 


( exp(iqn) + r^ r \q) exp {—iqn) (n > N) 


\ 


t ( r ' ) (q) exp (iqn) 


{n < —N) 


(24) 


(25) 


for right-side incidence, where 0 < q < ir and 


E = 2k cos q 


(26) 


is the real energy. The spectral reflection and transmis- 
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FIG. 5. (color online) Numerically-computed behavior of (a) 
spectral transmittance |t^(<?)| 2 , (b) reflectance |r^(q)| 2 , and 
(c) phase of t^(t) (left side incidence) in a heterojunction 
lattice with disorder at sites —N < n < N with N = 10 for 
10 realizations of disorder. The energies W„ are taken from a 
uniform distribution in the range (—k,k). The gauge field is 
h = 2.5. The dashed curves [almost overlapped with the solid 
ones in panels (b) and (c)] show the behavior of the spectral 
functions in the absence of disorder. In (d) the behavior of the 
spectral transmittance (q )| 2 in the absence of the junction, 
i.e for h = 0, is also shown for comparison. 


sion coefficients t^ l ’ r \q), r^ l,r \q) for left and right inci¬ 
dent sides can be determined following the general pro¬ 
cedure outlined in the Appendix A. The spectral coef¬ 
ficients are functions of the gauge field h and their val¬ 
ues as h —► 0 are denoted by tg’ r \q), ,r \q)- In the 
absence of the gauge field (,h = 0), the spectral trans¬ 
mission does not depend on the incidence side, i.e. one 
has tg\q) = ty\q) = t 0 (q), whereas for li / 0 a non¬ 
reciprocal behavior arises, i.e. t (l \q) ^ t^ r \q). The 
spectral coefficients are strongly affected by the two in¬ 
terfaces and by the different dispersion relations of the 
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lattices in the regions with and without the applied gauge 
field. In particular, a major role is played by the complex 
bias energy S between Hermitian and non-Hermitian re¬ 
gions. The case S = 0 is a rather simple one to handle 
analytically, however regrettably in this case the dou¬ 
ble heterostructure fails to realize robust and disorder- 
insensitive transport. In fact, assuming S = 0 the spec¬ 
tral problem for the amplitudes ip n [Eq.(4)] in the h ^ 0 
case can be mapped into the one with h = 0 after the 
substitution 

( 1 n<—N 

tp n = ‘t’n x \ exp (hn + hN) —N < n < N (27) 

^ exp(2 hN) n> N, 

i.e. one has 

E(j*n = ^(4*n +1 T 4*71— l) T V n (f) n . (28) 

For S = 0 the following simple relations are thus readily 
obtained 

A l ’ r \q) = r^ l ’ r \q) 

t^(q) = exp(2WV)t 0 (q) (29) 

&\q) = exp(—2 hN)t 0 (q). 

This means that, while the gauge field h makes the trans¬ 
mission non-reciprocal, amplifying forward-propagation 
waves and attenuating backward-propagating waves, it 
does not suppress the detrimental effect of disorder as 
compared to the Hermitian lattice. Such a result seems 
apparently to contradict the property of non-Hermitian 
transparency discussed in the previous section. The rea¬ 
son thereof stems from the fact that non-Hermitian trans¬ 
parency occurs for propagative waves (i.e. with a real 
Bloch wave number); however, for S = 0 the incident 
wave with propagative Bloch wave number q and energy 
E(q) = 2k cos(q) coming from the Hermitian lattice sec¬ 
tion excites, in the non-Hermitian interface region, two 
evanescent waves with complex wave numbers qi ,2 which 
are obtained from the relation (elastic scattering) 

2k cos q = 2kcos((7i i 2 — ih) (30) 

he. < 71,2 = ±q + ih. To excite propagative waves in the 
non-Hermitian interface region and thus to realize non- 
Hermitian transparency, we introduce a complex energy 
bias S, so that the condition (30) of elastic scattering 
changes as follows 

2k cos <7 = 2kcos(<7i,2 — ih) — iS. (31) 

For a given value of q, we may choose S so as <71 turns out 
to be real and positive, corresponding to the excitation 
of a forward-propagating wave in the non-Hermitian in¬ 
terface. We typically assume q = 7 r/ 2 , which corresponds 
to a propagative wave packet with largest group velocity 
and minimal distortion !32j . By imposing q = qi = 7t/ 2, 
from Eq.(31) one obtains 

(32) 


The site energies V±n at the interface sites n = —N and 
n = N are then chosen to obtain zero reflection and full 
transmission at q = 7t/ 2 in the absence of disorder, i.e. 
for W n = 0 (impedance matching). Their values can 
be readily calculated by imposing that Eq.(4) is satisfied 
by V’n = exp(— iqn) over the entire lattice at q = 7r/2. 
Assuming V n = — iS for —N < n < N and V n = 0 for 
\n\ > N one obtains 

V_n = —*«[ 1 — exp(— h)] , Vn = — iK[exp(/i) — 1]. (33) 

In the optical implementation of the Hamiltonian (2) 
shown in Fig.1(b), the imaginary energy bias S as well 
as V±n just correspond to the introduction of a suit¬ 
able level of optical loss (absorption or scattering) in the 
main microrings at — N < n < N. In practice, this can 
be achieved by a judicious doping of the dielectric mate¬ 
rial with atomic absorbers or by tailoring the scattering 
(radiation) losses (i.e. cavity Q factor) of the micror¬ 
ings by e.g. surface roughness control. The double het¬ 
erostructure lattice synthesized in this way realizes one¬ 
way transport which is robust against disorder in the 
non-Hermitian section of the junction. In fact, follow¬ 
ing a similar perturbative procedure as the one detailed 
in the Appendix B the asymptotic expressions for the 
spectral reflection and transmission coefficients, for both 
left and right incidence sides, can be derived in the limit 
h — 00 , which read explicitly 

t (Z )(g) = exp[2fqiV — iNn] —-|- O(e) (34) 

r«( 9 ) = -«p(2„JV)^=^+0<«) (35) 

t W («) = 0(e) (36) 

r M(q) = l + 0( e ) (37) 

where e = exp(— h). Note that at leading order in e the 
spectral coefficients do not depend on the disorder W n . 
Figure 4(a) shows the numerically-computed behavior of 
the spectral transmittance and reflectance for increasing 
values of the gauge field h in the absence of disorder 
(W n = 0), together with their asymptotic expressions 
given by Eqs. (34-37). The figures clearly shows that the 
heterostructure is highly non-reciprocal. In particular, 
at around q = n/2 transparency for forward-propagating 
waves, and full rejection of backward waves, is observed. 
It is also worth commenting on the behavior of the 
phase ipt of the spectral transmission t4 l \q) for forward- 
propagating waves, which is depicted in Fig.4(b). Note 
that, as h is increased, the group delay r = — ( d4>t/dq ) 
at around q = 7t/2 is non-vanishing and negative, indi¬ 
cating that a wave packet will travel the non-Hermitian 
region faster as h increases. Such a behavior can be sim¬ 
ply explained after observing that the group velocity v g 
of the wave packet in the non-Hermitian region of the 
junction is v g = 2k cosh h sin q, which is coshh times 
larger than the group velocity in the Hermitian regions 
[see also Fig.6(c) discussed below]. The impact of disor¬ 
der on transport along the junction is illustrated in Figs.5 


S = 2Ksinh(/i). 




and 6. Figure 5 shows the numerically-computed spectral 
functions for 10 different realizations of disorder in the 
same heterostructure of Fig.4 for h = 2.5. Note that the 
impact of disorder is fully negligible; such a result should 
be compared with the behavior of the spectral transmit¬ 
tance in the absence of the gauge field, i.e. for h = 0, 
which is depicted in Fig.5(d). The one-way transparency 
of the heterojunction for a Gaussian wave packet that 
crosses the disordered region is shown in Fig.6. While 
for h = 0 transmission is prevented for both left- and 
right-side incidence [Fig.6(a)], one-way transmission is 
observed for h ^ 0 [Figs.6(b) and (c)]. In particular, as 
h is increased, the wave packet travels faster in the non- 
Hermitian region [compare Fig.6(b) and (c)], according 
to the theoretical analysis. 


IV. CONCLUSIONS 

Edge states with topological protection that can prop¬ 
agate at the boundaries immune of disorder and scatter¬ 
ing are known to arise for electrons or bosons in two- 
dimensional systems subjected to magnetic or synthetic 
gauge fields. On the other hand, the effect of a magnetic 
field in one-dimensional systems, such as in a linear tight- 
binding lattice, is generally trivial and for such a reason 
it has received few attention so far. Long-range inter¬ 
action in one-dimensional systems can effectively mimic 
two-dimensional lattices, leading to fractal energy spec¬ 
tra, flat bands with localized edge states, and topological 
many-body states in the presence of synthetic gauge fields 
S3- However, the intrinsic one-dimensional nature of 
transport is not immune of Anderson localization. In this 
work we have presented a route toward the realization 
of robust one-way transport in one-dimensional lattices, 
based on the application of an imaginary gauge field. As 
shown in a pioneering work by Hatano and Nelson [21] . in 
a one-dimensional lattice with disorder application of an 
’imaginary’ (rather than real) magnetic field can prevent 
Anderson localization, with the appearance of a mobil¬ 
ity interval at the center of the band. In our work we 
have presented a detailed analytical study of the scatter¬ 
ing properties of one-dimensional lattices subjected to an 
imaginary magnetic field, and have highlighted the phe¬ 
nomenon of one-way non-Hermitian transparency: be¬ 
sides to prevent Anderson localization and restoring mo¬ 
bility, we have shown that for a sufficiently high imagi¬ 
nary magnetic field the structural imperfections and dis¬ 
order in the lattice become almost one-way transparent , 
i.e. they do not scatter back waves and do not substan¬ 
tially alter the spectral transmission of the lattice for 
one direction of propagation. Such a result stems from 
the fact that back-scattered waves are evanescent (rather 
than propagative) waves when an imaginary gauge field 
is applied. Interestingly, the spectral transmission of 
the non-Hermitian lattice can be mapped into the one 
of the corresponding Hermitian lattice, i.e. without the 
imaginary gauge field, but computed in the complex en¬ 


ergy plane. In particular, owing to complexification of 
the energy at large values of the gauge field the spec¬ 
tral transmittance becomes equal to one, even in the 
presence of disorder or lattice imperfections, thus real¬ 
izing non-Hermitian transparency. We have also shown 
that robust one-way transport and one-way transparency 
can be realized in a more realistic case of a double het¬ 
erostructure system, in which the magnetic field is ap¬ 
plied to a spatially-localized region of the lattice. Our 
results highlight important and rather unique effects of 
non-Hermitian transport and are expected to stimulate 
further theoretical and experimental investigations. In 
particular, the proposal to implement imaginary gauge 
fields in optical systems [32] [33], for example using a 
chain of coupled microrings [32] , would enable the obser¬ 
vation of important phenomena like non-Hermitian de- 
localization transitions and non-Hermitian transparency 
in heterojunction lattices. 

Appendix A: Spectral transmission and reflection 
coefficients 

In this Appendix we derive explicit forms for spectral 
transmission and reflection coefficients, corresponding to 
either left or right incidence sides, in a rather generic 
lattice described by the Hamiltonian T~L given by Eq.(l) 
in the text. The method is based on a rather standard 
transfer matrix approach. We assume that the lattice is 
asymptotically homogeneous and Hermitian, i.e. we will 
assume 

V n = 0 n< -N, n > N (Al) 

K n ,n+ 1 = Kn+I,n = « n < -N, n> N. (A2) 

A scattered state solution to Eq.(4) corresponding to left¬ 
side incidence has the asymptotic behavior 

exp(— iqn) + r^ ( q ) exp (iqn) (n < —N ), . 
t® (q) exp(— iqn) (n > N) 

where r^(q) are the spectral transmission and 

reflection coefficients for left-side incidence, respectively, 
and 0 < g < 7r/2 is the Bloch wave number of asymp¬ 
totic plane waves. Similarly, a scattered state solution 
to Eq.(4) corresponding to right-side incidence has the 
asymptotic behavior 

I f &\q)exp(iqn) (n<-N),., 

n \ exp (iqn) + r^(q) exp {—iqn) {n > N) 

where U r ) (g), r ^ (q) are the spectral transmission and re¬ 
flection coefficients for right-side incidence, respectively. 
The difference equation (4) can be cast in the following 
matrix form 

(V )-"-(&) (A5) 

and 

(A6) 
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lattice site n lattice site n lattice site n 


FIG. 6. (color online) Propagation of a Gaussian wave packet (snapshots of |c n (t)| 2 / l c ™(/)| 2 ) i 11 a heterojunction lattice 

with disorder for (a) h = 0 (Hermitian limit), (b) h = 0.5, and (c) h = 2. The energies W n at sites —N < n < N (with N = 8) 
are taken from a uniform distribution in the range (— k, k) (see the inset). Upper panels refer to left-side incidence (carrier 
Bloch wave number of the incident wave packet q = n/2), whereas lower panels refer to right-side incidence (carrier Bloch wave 
number of the incident wave packet q = —tv/2). 


where we have set 


_ ( {E Vn )/ttn,n+l ^n,n—l/ttn,n+l \ 

y 1 0 J 


(A7) 


Note that M n = M n in the absence of the gauge field, 
i.e. for h = 0. 

From Eqs.(A3), (A4) and (A9) it then follows 


(A Vn) / Kn,n— 1 ^n,n+l/^n,n—1 

1 0 

and E = 2k cos q. Hence one has 

(V )-“(&,) (A9 » 

and 

(^r)= p U«) (Aio) 

where we have set 

Q = Mn x Mn~i x •••• x M-n+i x M-n- (All) 
and 

V = M-n x M-n+i x .... x Mn—i x Mn- (A12) 




( (q) exp[-iq{N + 1)] \ _ f Qu Qu ) 

\ (q) exp(-iqN) ) ^ Q 21 Q 22 J 

( exp(iqN) + (q) exp(—iqN) VaiA 

\ exp[*g(lV + 1)] + (q) exp[-iq(N + 1)] J 


exp[iq(N + 1)] + r^ (q) exp [—iq(N + 1)] ^ _ 
exp(iqN) + r^ ( q) exp(-iqN) 


( Qu Oil t < - r '>(q)exp(-iqN) 

\ Q 21 Q 22 J \ t( r \q) exp[— iq(N + 1)] J 
where Qik are the elements of the 2x2 matrix Q. Equa¬ 
tions (A13) and (A14) can be solved for A l \q) and 
t {r \q), yielding 


(q) = exp[iq(2N + 1)] 
t^ r \q)=exp[iq(2N+l)} 


Q 21 exp(-iq) - Q 12 exp(iq) + Q 22 - 
Q 11 exp (iq) - 0 22 exp(-iq) + Q 12 - 
2 i sin q 

On exp (iq) - Q 22 exp (-iq) + 0r 2 - 


On 

Q 21 

(A15) 

Q 21 

(A16) 


Similarly, from Eqs.(A3), (A4) and (A10) one has 


f exp [iq(N + 1)] + r^(q) exp [—iq(N + 1)] ' 
\ exp(iqN) + r^{q) exp(— iqN) 

( Vu V l2 \ ( t( l \q) exp(-iqN) 

\ V 21 V 22 J y t m (■ q ) exp[— iq(N + 1)] 


(A17) 


( t^ r \q) exp[—iq(N + 1)] \ = ( Vn Vi 2 \ 

\ t^ r \q) exp(—iqN) J \V 2 i V 22 J 

( exp(iqN) + {q) exp(-iqN) \ . . 

\ exp [iq(N + 1)] + A r \q) exp [~iq{N + 1)] J 
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which can be solved for r^ r \q), t^(q), yielding 


r^ r \q) = exp[iq(2N + 1)] 
t®(q) = exp[iq(2N+l)] 


V 2 i exp (-iq) - V l2 exp (iq) + V 22 
Vn exp (iq) - V 22 exp (-iq) + V\ 2 
2 i sin q 

Vn exp(iq) - V 22 exp (-iq) + V\ 2 


Vn 
- V 21 

(A19) 

v 21 

(A20) 


Appendix B: Asymptotic behavior of spectral 
coefficients 

In this Appendix we prove the asymptotic form of 
the spectral coefficients t(q) = to(q — ih ) and r(q) = 
exp(—2 Nh)r 0 (q — ih) given by Eq.(20) in the text. To 
this aim, let us notice that the expressions of t$(p) and 
ro(p) are given by Eqs.(9) and (10), where p = q — ih 
and the matrices Q, V are defined by Eqs.(11-13) given 
in the text. After setting e = exp (—ip) = exp (—iq — h), 
we consider the asymptotic forms of Q and V as h —► 00 , 
i.e. e — > 0, assuming V n /n of order ~ 1 (or smaller). 
Since V is obtained from the expression of Q after the 
substitution V n —► V_ n , we can limit to compute the 
asymptotic form of Q. To this aim, we note that one 
has Q = e -( 2Ar + 1 ) 7 z with TZ of order ~ 1. Owing to the 
explicit dependence of to(p), ro (p) on e = exp (—ip) [see 
Eqs.(9) and (10)], it is enough to compute the asymptotic 
expressions of the coefficients 7 Zik of TZ up to the order 


~ e for (i,j) ^ (1,2), whereas the asymptotic expression 
of 1Z\ 2 should be pushed up to the order ~ e 2 . At such 
an order of approximation, after some calculations one 
obtains 


Q 


1 

e 2N+l 


1-eE 


N 

1——N 

€ 


Yl 

K, 


e + e 


: sr^N-l Vi_ 

2^i=-n k 

0 


( B1 ) 

The asymptotic form of the matrix V is obtained from 
Eq.(Bl) after the change V n —> V- n . Substitution of 
Eq.(Bl) (and a similar expression for V) into Eqs.(9) 
and (10) yields 


t(q) = t 0 {q - ih) = 1 + O(e) 

(B2) 

r(q) = r 0 (q — ih) exp(—2 hN) 


_ eV/v ex p( 2 iqN) + 0(e 2 ) 

(B3) 


which prove Eq.(20) given in the text. 
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